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Vector Area Theorem mapping in crystals and polarization stability of 

SIT-solitons 
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The stability of polarization, areas, and number of self-induced transparency (SIT)-solitons on an output 
from the LaFz ■ Pr 3+ crystal is theoretically researched versus of the polarization direction and the area of 
the input linearly polarized laser pulse. For this purpose the Vector Area Theorem is rederived and two- 
dimensional Vector Area Theorem map is obtained. The map is governed by the crystal symmetry and take 
into account directions of the dipole matrix element vectors of the different site subgroups of optically excited 
ions. The Vector Area Theorem mapping of the time evolution of the laser pulse allows to highlight soliton 
polarization properties. 
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For an isotropic medium stability properties of SIT- 
solitons are determined by the Area Theorem. The 
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Area Theorem is the name given to a theoretical re- 
sult that governs the coherent nonlinear transmission of 
short light pulses through isotropic materials that have 
an absorption resonance very near the frequency of the 
frequency of incident light, effectively two- level media. 
In 1967 McCall and Hahn [1] identified a new param- 
eter (called "Area" and denoted by 9) of optical pulses 
travelling in such media, and then predicted that obeys 
the simple equation 
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.where a is the attenuation coefficient for the material. 
^ (— | The two most striking consequences of the Area Theo- 
Cl^rem are: (i) pulses with special values of Area, namely 
-> all integer multiples of n, are predicted to maintain 
the same Area during propagation, and (ii) pulses with 
other values of Area must change in propagation until 
their Area reaches one of the special values. This prop- 
erty can be shown to be unstable for the odd multiples, 
but the even multiples enjoy the full immunity of the 
theorem. The Area Theorem (1) was derived for an 
isotropic material in which the dipole matrix element 
vector of any ion is parallel to the electrical field vec- 
tor of the light pulse. By contrast, the direction of the 
dipole matrix element vector of any Pr 3+ ion in LaF^ 
is not dependent from the electrical field vector of the 
light pulse. So it is necessary to rederive Area Theorem 
taking into account directions of the dipole matrix el- 
ement vectors of the different subgroups of Pr 3+ ions. 
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The Pr 3+ ions in a LaF^ unit cell can replace La 3+ 
in six different types of sites (±a,±/3, ±7). The local 
environment of any of them has ^-symmetry. The six 
local CVsymmetry axes are located in the plane normal 
to the C3 axis and make the angle of 2ir/6 in this plane 
(Fig. 1). The electrical dipole matrix element vector of 
the Pr 3+ ion (optical transition Ti — ► Tx) 

•,6 (2) 
-symmetry axis 

e z ■ Bj = 0. (3) 



Pj =P e jr--3 = Ir- 
is directed [2] along of the local C2 

e j = ( c os(jy),sin(jy)), 



where ej is unit vector along the C2j axis. Here axis Z 
is directed along the C3-axis and axis X along a-axis. 
We define the Vector Area of light pulse as 
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where E(z, t) is vector amplitude of a light pulse, p is 
electrical dipole matrix element, h is the Plank's con- 
stant. Taking into account (2) and (3) and using ar- 
guments as Lamb [3] we can write the Vector Area 
Theorem as follows: 
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if a light pulse propagation is along C3-axis. 
the linear attenuation coefficient for LaF^ : 
can be seen from (5), as © — > , that 



d® 

dz 



(5) 

Here a is 
Pr 3+ . It 
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as it must for a small pulse Area. 

Equating the right part of the equation (5) to zero 
we can find special values of the Vector Area where 
d@/dz = 0. But it to make much more obviously and 
more easy from the graphical representation. We can 
rewrite eq (5) as 
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And the problem is reduced to a determination of 
points in a two-dimensional plane, in which the function 
Y^cos(0 -ej) has extremes. The circles and triangles in 
Fig. 2 give the contour plot of this function. We easy 
find three types of special values of the Vector Area, 
namely 



© c = m©+ + n® 



® u = ® c + u, 



S = C + Sj, 



(8) 



(9) 



(10) 



are predicted to maintain the same Vector Area during 
propagation. Here m and n are arbitrary integers and 
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where the unit vectors e 3 (3) and kj are directed along 
and between the C2-axes accordingly: 
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k 3 = ( cos (i-g- - g ). sin (i-g- - g)) 
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These peculiar points (8,9,10) give rise to a two- 
dimensional lattice in a ©-phase plane with basis vec- 
tors © + and (11) as you can see in Fig. 2. The 
unit cell of the lattice is determined by symmetry of the 
crystal. It is a regular hexagon. The hexagon's centers 
arc © c (8) (centers of the circles in Fig. 2). As mea- 
sured from the hexagon center, coordinates of the six 
tops of the hexagon are Uj (13) (centers of the trian- 
gles in Fig. 2), coordinates of middles of the sides of a 
hexagon are Sj (12). As is easy to see from Eqs. (5,7) 
and definitions (8,9,10), that in a neighborhood of these 
singular points the Vector Area behaves as 
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if — ® s is directed along the side of a hexagon, and 

(18) 
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if — S is directed perpendicularly to the side of 
a hexagon. Therefore, for an absorbing (amplifying) 
medium with a > (a < 0), the points (8) are type 
of a stable (unstable) knot, the points (9) are type of 
an unstable (stable) knot and the points (10) are type 
of a saddle in the 0-phase plane. In further we shall 
explore a case of the absorbing medium with a > 0. 
If an input pulse Vector Area falls inside an unit cell 
then the Vector Area must change in propagation until 
it reaches the unit cell center. If an input Vector Area is 
not equal to (9-10) and falls on a side of a hexagon then 
the Vector Area must change in a propagation until it 
reaches the middle of the gexagon side. It is necessary 
to mark that the Vector Area Theorem map (Fig. 2) 
allows easily to predict only the sum of the pulse vector 
areas on an output from the sample. To determine the 
number of the output SIT-solitons and their polariza- 
tions and areas we should solve the system of coupled 
Maxwcll-Bloch equations. 

The input pulse, which Vector Area is directed be- 
tween the crystallographic axes and is equal, for exam- 
ple, to + , excite only four (±a, ±7 ) ion subgroups. 
It is 27r-pulse for these ions. This pulse does not excite 
±/3-ions, because 0+_Le± ) g. To the similarly previous 
case, the input pulse, which Vector Area is equal to 
, is 27r-pulse for (±a, ±/3 )-ion subgroups and does 
not excite ±7-ions. If the input Vector Area is parallel 
to the + (0-) and fall inside the unit cell C = m® + 
(® c = m0_ ), then the time evolution of pulse, as it is 
easy to show, may be described by the inverse scattering 
method. The input pulse is splitting up on output, as for 
an isotropic medium, into m SIT- solitons, which Vector 
Area of any of them is 0+ (©-)■ In further we shall 
refer to these solitons as + -solitons and 0_-solitons. 
If input Vector Area is not parallel to the + (©-) 
but fall inside the unit cell ® c = m® + (0 C = m©_ ), 
then, as show numerical calculations, input pulse also is 
splitting up into m © + (@_)-solitons on output. 

If an input pulse Vector Area is directed along a crys- 
tallographic axis, for example the axis a, and is equal 
to 



©o = 0+ + 0- 



(19) 



Vector Area Theorem mapping 
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then all (±a, ±/3, ±7) ions are excited. The input pulse 3. G. L. Lamb, in: Elements of Soliton Theory, John Wil- 
is 27r-pulse for (±/3, ±7 ) ion subgroups and 4p-pulsc ley and Sons, New York, 1980. 
for (zba ) ions. The time evolution of the pulse is not 
described by the inverse scattering method. The nu- 
merical calculations have shown that if an input Vector 
Area 0;„ is parallel to ©o and fall inside the unit cell 
C = m&o, then the input pulse is splitting up into 
m SIT-solitons, which a Vector Area of any of them is 
©o • In further we shall refer to these solitons as ©o- 
solitons. Let an input Vector Area is not parallel to 
@o and fall inside the unit cell C = to© . Then, as 
you can see in Figs. 2-3, a small change of an input 
pulse polarization reduces to that each of the ©o-soli- 
tons is splitting up into + -and @_-solitons. Therefore 
a number of solitons and its polarization strongly de- 
pend on a direction of a vector @j„ concerning a crys- 
tallographic axis. This deduction also is valid and in 
a generally case when the input Vector Area full down 
inside the unit cell C = to© + + n® , where m =/= n. 
It is so because the unit cell center coordinates may be 
rewriting as C = (to — n)® + + n®o if m > n or as 
C = (n — m)0_ + m&o if n > m. At first there are 
(to — n)0 + -solitons for m > n, or (n — m)0_-solitons if 
n > m on an output. After that the number of solitons 
appearing on an output depends on a direction of the 
vector 0i„ — (to — n)0 + or 0j n — (n — to)0_ concerning 
a crystallographic axis. In a stable case the output soli- 
tons are + -solitons and 0_-solitons and its number is 
(to + n) . 

For an amplifying medium + -solitons and 0_ -soli- 
tons are unstable so the polarization of the output soli- 
tons must be directed along crystallographic axis in a 
stable case. 

It is necessary to mark that for circular polarization 
of laser pulse the Area Theorem is (1) as in a case of 
isotropic medium. 

To summarize, we have shown on an example of the 
model system LaF 3 : Pr 3+ , that the Vector Area map- 
ping of the pulse time evolution during a propagation is 
effective method to analyze the polarization properties 
of solitons. 
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00-02-16510. 
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Fig. 1. The directions of the local ^-symmetry 
axes for the different Pr 3+ ion sites in the plane normal 
to the C3- axis of the LaF^ : Pr 3+ crystal. 

Fig. 2. Vector Area Theorem map. The ®/2ir pro- 
jections to axes X and Y are plotted on axes X and ac- 
cordingly Y. The vector si (12) and Ui (13) are shown 
in an upper of the figure. There are the basis vectors 
©_l_ and (11) and the unit vectors along of the local 
C2-symmetry axes (+a, —/3, —7) in a down. Vector co- 
ordinates of some unit cell centers are also shown. Bold 
lines and 1 are mappings of the time evolution of laser 
pulses which the input Vector Area is mod(© ) = 47r 
and the angles between the directions of the Vector Area 
and the crystallographic axis a are and —1 degrees ac- 
cordingly. In this case aL = 20, where L is the sample 
length and a is the attenuation coefficient. The bold 
line 2 is the mapping of the time evolution of the laser 
pulse which the input Vector Area is mod(2& ) = 8n 
and the angle between the directions of Vector Area and 
the crystallographic axis a is +1 degree and aL = 40. 
The bold line 3 is as 2 but aL = 80. 

Fig. 3. The time evolution of the amplitude of the 
laser pulses on the output of the sample. Values of the 
input Area, the angles between the directions of the 
Vector Area and the crystallographic axes a and the 
parameter aL for the curves 0,1 and 2 are the same as 
for the curves 0,1 and 2 in Fig. 2. The parameter value 
aL for the degree curve in down of the figure is as for 
curve 2. The dot line is input pulse, t p is input pulse 
duration. The Vector Area Theorem mapping (curves 1 
and 2 in Fig. 2) allows easily to spot the polarizations 
and the areas of the solitons in curves 1 and 2 in this 
figure. 




X-proection of Vector Area (2ti) 
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input Area = 4jt 





